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We study the Bethe approximation for a system of long rigid rods of fixed length k, with only 
excluded volume interaction. For large enough k, this system undergoes an isotropic-nematic phase 
transition as a function of density of the rods. The Bethe lattice, which is conventionally used to 
derive the self-consistent equations in the Bethe approximation, is not suitable for studying the 
hard-rods system, as it does not allow a dense packing of rods. We define a new lattice, called 
the random locally tree-like layered lattice, which allows a dense packing of rods, and for which 
the approximation is exact. We find that for a 4-coordinated lattice, fc-mers with k > 4 undergo 
a continuous phase transition. For even coordination number q > 6, the transition exists only for 
k > k m i„(q), and is first order. 

PACS numbers: 64.70.mf,05.50.+q,64.60.Cn 



I. INTRODUCTION 

The study of the ordering transition in systems with 
only excluded volume interactions has a long history. De- 
pending on the shape of the molecules involved, these 
systems can exhibit different ordered states and undergo 
transitions between them. Many different shapes have 
been studied in literature: e.g., hard spheres hard 
squares 0,0, hexagons trianglesjl , tetrominoes @, 
rods Q, banana-shaped molecules [Jetc. A system of 
thin rigid cylindrical molecules in solution was shown by 
Onsager to undergo a first order phase transition from 
an isotropic to a nematic phase at sufficiently high den- 
sities, in the limit of large aspect ratio 0. Flory studied 
the hard rods problem on a lattice , still allowing the 
rods to have continuous orientations, and using a mean- 
field approximation again found a isotropic-nematic tran- 
sition. Zwanzig studied a system of hard rods in the con- 
tinuum, but restricting their orientations to a discrete set 
[llj , also finding a discontinuous transition to a nematic 
phase as the density of rods is increased. The relation 
between continuous and discrete models is discussed in 
Ref. [13 . 

For continuous models, in dimensions d > 3, if the as- 
pect ratio of the rods is high enough, the existence of an 
isotropic-nematic transition is well accepted Q. In two 
dimensions, in a system of hard needles in a plane, the 
Mermin- Wagner theorem [l3j disallows long range orien- 
tational order, but at sufficiently high densities, the sys- 
tem shows a Kosterlitz-Thouless type ordering transition, 
and the orientational correlations decay with distance as 
a power law .14 , UH . If the rigid rods are placed on a lat- 
tice, then only a discrete set of orientations is possible, 



and the existence of a long-range ordered phase is not 
disallowed by the Mermin- Wagner theorem. Rods occu- 
pying k consecutive sites along any one lattice direction 
will be called fc-mers. For dimers (fc = 2), it can be shown 
rigorously that orientational correlations decay exponen- 
tially, except in the limit when all sites are occupied by 
rods. In the case of full coverage, the orientational cor- 
relations have a power law tail in all dimensions d > 2 

USE!- 

The question whether the lattice model of rigid rods 
with k 3> 2 undergoes a phase transition remained unset- 
tled for a long time [ll|. Clearly, the number of ways of 
covering the planar square lattice fully with fc-mers grows 
exponentially with the size of the system, and most of 
these configurations are not ordered. Thus, at the maxi- 
mum packing density, which was presumed to favour or- 
dering most, the system shows no ordering. Recently, 
Ghosh and Dhar argued that fc-mers on a square lattice, 
for fc > 7, would undergo two phase transitions, and the 
nematic phase would exist for only an intermediate range 
of densities p* < p < p\, where p? ~ 1/fc [l]|- Subse- 
quently, Fernandez et. al. |20l - |23j . through extensive 
simulations and approximate methods, showed that the 
square lattice transition is Ising like, or equivalently in 
the liquid-gas universality class 24|. On the triangular 
and hexagonal lattices, there are three degenerate fully 
aligned ground states, and the transition was shown to 
be in the q = 3 Potts model universality class [13, HH ■ 
The existence of a transition to a state with orientational 
order was shown rigorously for a polydispersed hard rod 
system on a two dimensional square lattice (25| . The sec- 
ond transition from nematic to unordered state at p* 2 is 
much less understood. In Ref. [19(, a variational estimate 
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of the entropy of the nematic and ordered states suggests 
that 1 — p* 2 should vary as 1/fc 2 for large k, but not much 
is known about the nature of transition. Linares et. al. 
estimated that 0.87 < p* 2 < 0.93 for k — 7, and pro- 
posed an approximate functional form for the entropy as 
a function of the density (23|. 

Given the paucity of exact results about the second 
transition, it seems worthwhile to investigate this prob- 
lem using well-known approximation methods, like the 
Bethe approximation (BA). However, we realized that 
there are problems in the application of the conventional 
BA techniques, even to the first isotropic-nematic tran- 
sition. In particular, packing fc-mers on the Bethe lattice 
with a density close to one seems difficult, and the na- 
ture of the high-density phase on the Bethe lattice is not 
very clear. Another way to study BA is to work with 
random graphs of fixed-coordination number [26-291] . In 
this case, there is no surface, but in this case as well, one 
does not get full coverage in the limit when the activity 
of fc-mers tends to infinity. This has led us to introduce 
a new lattice, to be called the random locally tree-like 
layered lattice (RLTL lattice). We study /c-mers on the 
RLTL and show that there is no difficulty in defining 
the high-density phase for the A:-mers and that the BA 
is exact on this lattice. In fact, we find many ways to 
cover the lattice with k-mers. There is a finite entropy 
per site in the fully packed limit, and in a typical fully 
packed configuration, all k-mers are not all in the same 
orientation. 

The rest of the paper is organised as follows. We begin 
with a recapitulation of the conventional approaches to 
the Bethe approximation, and the Bethe lattice in SecUTl 
and derive self-consistent equations between different cor- 
relation functions. In Sec. IIII| we describe the RLTL 
lattice where every site in the lattice has the same even 
coordination number q. In Sec. IIV1 we derive the exact 
annealed partition function of the problem on this lattice, 
and show that the resulting self-consistent equations are 
the same as obtained earlier for the Bethe lattice. In 
Sec. El we analyse the behaviour of these equations for 
the RLTL with coordination number 4, and show that it 
undergoes a continuous transition for k > 4. In Sec. IVI| 
we discuss the case q > 6. We find that in this case, the 
system undergoes a first order transition for k > k m i n {q). 
Section [VIII contains some concluding remarks. 



II. fc-MERS ON A 4-COORDINATED BETHE 
LATTICE 

The BA was initially proposed as an ad-hoc decoupling 
approximation (the pair approximation) [30L [3l| . It was 
realized that the approximation becomes exact for the 
Bethe lattice, which is the part of a uniformly branch- 
ing tree graph far away from the surface 32|, 33 1 . This 
is important, as this ensures that the approximation can 
not give rise to unphysical consequences, like non-positive 
specific heat. Also, there is a consistent scheme, and two 



different answers for same quantity cannot be obtained. 
Given that in a uniformly branching tree, most of the 
sites are a finite distance from the surface, disentangling 
the surface contribution from that of sites deep inside re- 
quires some care [33 - l36| . There are some prescriptions 
for separatin g th e surface contribution from bulk given 
in literature [37|. A rather careful and detailed discus- 
sion for the case of the Ising model on the Bethe lattice 
is given in Chapter 4 of Ref. [38| . One considers corre- 
lation functions deep inside a Cayley tree, whose values 
are independent of the boundary conditions on the sur- 
face of the tree. This can be realized for all non-zero 
external magnetic fields. Then thermodynamical quanti- 
ties like the free energy are obtained by integrating these 
correlation functions with respect to suitable conjugate 
field. 

A more recent approach to BA is to treat it as the free 
energy for a random graph on N sites, where each site 
in the graph has the same degree. One can show that in 
the limit N tends to infinity, for a randomly picked site 
in the graph, the length of the shortest loop increases as 
\ogN . Then, the graph looks like a loop-less of uniform 
coordination number, up to a distance of order log N. In 
the limit of large N, it looks like a Bethe lattice, except 
that there is no surface. This makes study of the different 
models on the Bethe lattice numerically feasible [2H4291 ] . 

We first derive the self consistent equations satisfied by 
densities for the model of fc-mers on a Bethe lattice of co- 
ordination number four in the conventional way. Gener- 
alisation to higher even coordination number is straight- 
forward. A Bethe lattice corresponds to the core of the 
Cayley tree with the same coordination number. The 
sites of the Cayley tree may be ordered by its genera- 
tion, starting with m — for the sites on the surface 
and ending at m = M for the central site of a tree with 
M generations. Each site of the Cayley tree, other than 
those on the surface, has 4 bonds attached to it. Two 
of them will be said to be of type X and the remaining 
two to be of type Y. A fc-mer occupies k — 1 consecutive 
bonds ( and k sites) of the same type. We associate an 
activity Z\ with an X-type fc-mer (referred to as X-mer 
in the following) and an activity zi with a Y-type £;-mer 
(Y-mer). No two A;-mers may overlap, so that each site 
of the lattice has at most one &;-mer passing through it. 
Fig- H shows a four-coordinated Cayley tree occupied by 
three 4-mers. 

The Cayley tree is composed of four rooted subtrees 
which are connected to the central site. We define partial 
partition functions for these rooted subtrees with fixed 
configuration of the root bond. We specify the type of 
the root bond and the configuration of the root bond. 
The root bond may be empty, or if it is occupied by 
a &;-mer, then we specify the number of sites in earlier 
generations already occupied by this /c-mer. There are 
Ik such partial partition functions, which we denote by 
fjxj and g Vt j, where the first subscript denotes the type of 
the root bond and j = 1, . . . , k — 1 is the number of sites 
that are currently part of the fc-mer. j — corresponds 
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FIG. 1. A four-coordinated Cayley tree with four generations 
is shown. Two X-mers (horizontal thick solid lines) and one 
Y-mer (vertical thick solid line) with k = 4 are placed on the 
tree. The weight of the configuration is z\z2- 



to the case when the root bond is empty. 

The recursion relations obeyed by the partial partition 
functions are obtained by building a subtree with one 
additional generation by connecting 3 subtrees to a new 
root bond and new root site. This process is illustrated 



m Fig. [2j We attach a weight z x and z 2 with each 
occupied site of and X-mer and Y-mer respectively. The 
recursion relations are (see Fig. |3]) 
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where the prime denotes partial partition functions of 
subtrees with one additional generation of sites. 

The partial partition functions are multiplied by the 
appropriate activity each time a fc-mer grows such that 
the weight of a fc-mer is z\ or zi depending on its type. 
The partial partition functions grow exponentially with 
the number of iterations. To calculate the densities in 
the core of the tree (Bethe lattice), it is useful to define 
ratios of partial partition functions 



Rij — — i — x,y and j = 0, . . . , k — 1, 

9i,0 



(2) 



FIG. 2. Building a subtree with M + l = 4 generations of sites 
by connecting q ~ 1 = 3 subtrees with M = 3 generations of 
sites to a new root site and bond (denoted by dashed lines). 
If the root bond of the new subtree is in the Y-direction, then 
two of the existing subtrees will have their root bond in the 
X-direction. 
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FIG. 3. Diagrammatic representation of Eq. Q for g' y j, j 



such that R x ,o — Ry,o — 1- The recursion relations 
obeyed by the ratios are easily derived from Eq. ([T]): 
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D, 
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j = l,...,fc-l, (3a) 
i = l,...,fc-l, (3b) 



where 



A— 1 



D x = 1 + z/ R x ,k~i + z\l ^ /''./.;/''./ /. i ;• (4a) 

fc-1 

D a = 1 + z^ k R Vi k-i + z\^ k ^ RxjRx^-i-j- (4b) 
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We study the fixed points of these recursion relations. 
The fixed point corresponds to distance of the surface 
sites very large, or equivalently to sites in the deep inte- 
rior of the Cayley tree. Also, if only one single stable 
fixed point exists, the boundary conditions at surface 
sites (which determine the initial values of the recur- 
sions), have no effect on the behaviour at points deep 
inside (Hill. 

A natural choice of the starting weights for the recur- 
sion equations are such that all fc-mers are completely 
contained within the Cayley tree. This corresponds to 



At the fixed points, the densities simplify to 



1 Ik 

Jx,o = g y ,o = 1, 9x,i = Z\ , and g y>1 — ^ 2 
9x.j = 9y,j = for all j > 1. However, for these choices 
of initial conditions, the recursion relations under iter- 
ation converge to a fixed point that corresponds to the 
isotropic system because the initial conditions themselves 
are symmetric. The isotropic fixed point is stable with 
respect to perturbations where the initial conditions arc 
still symmetric with respect to the two directions. How- 
ever, it is unstable to asymmetric perturbations. 

For generic initial conditions, we find that on repeat- 
edly iterating the recursion relations, the R's converge 
to a stable fixed point R* , which may be determined by 
iterating Eq. fl3J. Then, 
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a? x , j = 0,. . .,k - 1, (5a) 



R* 



= a yl j = 0,...,fc-l, (5b) 



where the variables a Xl a y satisfy the equations 
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Knowing the fixed point solution, the density at the 
central site may be calculated. The grand canonical par- 
tition function of the system on the Cayley tree is given 
by 

fc-i 

— 2 2, l/k 2 \ ^ 

- = 9x.o9y.O + Z l 9y,0 2^9x,j9x,k-l-j 

3=0 

fc-1 

+ 4 9x,0 9y,j9y,k-i-j- (7) 

3=0 

Then, the densities p x (p y ) of sites that are part of X- 
mers (Y-mers) are given by 
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Eliminating a x , a y from Eqs. @ and Eqs. we 
obtain 



Zl(l - Px - Pyf 
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fe-1 
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Z 2 (l- Px -Py) k =^(l-- 1 ^P y 
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The Bethe lattice solution described above is not very 
satisfactory. In particular, in the limit of large z\ and zi 1 
one gets the fraction of sites occupied by A:-mers tends to 
1. However, it is easy to see that this can only be achieved 
with a very special choice of boundary conditions at the 
surface of the tree. 

Also, if we consider a tree with coordination number 
q > 6, in the large activity limit, then at any root vertex 
of a subtree, the process of exchanging different branches 
is clearly a symmetry operation on the lattice. It is easy 
to see that the presence of this local symmetry implies 
that there can be no nematic order in the deep inside 
region, with one type of bonds preferentially occupied. 
This makes this approach unsuitable for the studies of 
the hard rod problem. 



III. THE RLTL LATTICE CONSTRUCTION 

For simplicity, we discuss the random lattice with co- 
ordination number q = 4. Generalisation to other co- 
ordination numbers is straight forward. We consider a 
set of M layers, numbered from 1 to M, with N sites in 
each layer. A layer m is connected to the adjacent layer 
rn — 1 by N bonds of type X and bonds of type Y. 
The connections are made by randomly pairing each site 
of the mth layers with exactly one site in the (m — l)-th 
layer with an X-bond, similarly randomly pairing each 
site with a site in the neighbouring layer using Y-bonds. 
The total number of such possible pairings is (iV!) 2 . This 
is illustrated in Fig. 

We can impose open or periodic boundary conditions. 
For a 2q coordinated lattice, there are (N\) qM different 
possible graphs, in the case of periodic boundary condi- 
tions (layer M is connected to layer 1), and (AT!) 9 ^ -1 ^ 
different possible graphs for the case of open boundary 
conditions. We consider annealed models on this lattice, 
thus we average the partition function over all possible 
configurations of the bonds. 

We can associate different degrees of freedom with the 
vertices, and consider statistical mechanical models on 
these lattices. For example, we can attach an Ising spin 
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The correlation functions can be denned as usual. Con- 
sider the example of the Ising model defined above. We 
consider two sites i and j, and consider the two-point 
correlation function (SiSj). Since this correlation func- 
tion has to be averaged over all assignments of bonds, it 
can depend only on the difference in layer numbers of the 
sites. In particular, it has the same value for all sites i 
and j in the same layer. Thus, the expected factorisation 
property of the correlation functions, which is the essence 
of the Bethe approximation, is built into the definition 
of the lattice. 



FIG. 4. Schematic diagram of the random lattice with N = 6 
sites per layer and coordination number q — 4. A possible 
configuration of bonds between layers m — 1 and m is shown, 
with the solid and dotted lines being bonds of two different 
types. 



Sj to each vertex j of the lattice, and define the Hamil- 
tonian to be H = —JJ2 nn &iSj, where the sum is over 
all nearest neighbour pairs. For a particular realization 
1Z of bonds, let the partition function be Z-ji(M, N). We 
average the partition function over all possible config- 
urations of bonds. As we are averaging the partition 
function, and not the logarithm of the partition func- 
tion, this is similar to the annealed average over different 
bond configurations. Thus, 



f 



Z av (M,N) = — J2Zn(M,N), 



(11) 
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where Nji is the number of different bond configurations 
on the lattice. 

To take the thermodynamic limit, we let M and N tend 
to infinity. The mean free energy per site / is defined as 
the limit 



/ 



-kT lim 

M,iV-K 



f 



MN 



log Z av (M,N). (12) 



We note that each site, say in layer m, has exactly 
one X-bond, and one Y-bond connecting it to sites in the 
layer m + 1. It can happen that these bonds fall on the 
same site. However, the fractional number of such cases is 
only 1/N , and thus, the expected number of loops of size 
2 on this lattice is M. For M even, the graph is bipartite, 
and there are no loops of odd perimeter. It is easy to 
verify that the expected number of loops of perimeter 4 
on this lattice is 5M. In general, the number of loops of 
perimeter I per site of the lattice varies as A /N, where 
A is the self-avoiding walk growth constant [39( . 

We, thus, see that for large N, there are very few short 
loops in the lattice. For a randomly picked site, the size 
of the shortest loop going through that site is of order 
log TV, and this goes to infinity, as N goes to infinity. 
Since there are very few short loops, the structure of the 
lattice locally is that of a regular branching tree, and it 
locally looks like the Bethe lattice. 



IV. fc-MERS ON THE RLTL LATTICE WITH 

9 = 4 

We consider a fc-mer model on the RLTL lattice. A k- 
mer occupies (k — 1) consecutive bonds of the same type. 
As earlier, we associate an activity Z\ with an X-mer and 
an activity Z2 with a Y-mer. 

Let x m be the number of X-mers with topmost (1 st ) 
site in the m-th layer, and y m is number of Y-mers with 
topmost site in the m-th layer. We will denote by X m 
and Y mi the number of sites in the roth layer, occupied 
by X-type, and Y-type bonds respectively, but where the 
site is not the topmost site of the fc-mer. Clearly, we have 



A',, 



fc-i 



fc-i 



Yn, 



Vm-j ■ 



(13) 



We can adopt the convention that x m = y m = 0, for 
m < 0. Then the above equation holds for all m, 1 < 
m < M. Note that if we want the fc-mers to be fully 
contained in the lattice, we must also have x m = y m = 0, 
for m > M - k + 2. 

To calculate the partition function, consider the oper- 
ation of adding an additional layer. We, thus, specify 
the full set of 2M values {x m ,y m } for all m. The total 
statistical weight of configurations that contribute to a 
particular set {x m ,y m } will be denoted by C({x m ,y m }). 

We calculate C({x m ,y m }) recursively. Let us imagine 
that we have constructed the configuration of the lattice 
to layer m—1, and now we add the layer m. We sum 
over the configurations of bonds, and that of X-mers and 
Y-mers on these bonds at the same time. 

(i) We note that in the (m — l)-th layer, there are X m 
sites that are occupied by X-mers which protrude to layer 
m. The first of these sites is connected to a randomly 
picked site in the lower layer in N ways using an X- 
bond. Then this site in the lower layer is occupied by 
the extension of that X-mer. The next can be connected 
to one of the unoccupied sites in the lower layers in only 
{N — 1) ways. Similarly the third, and so on. The number 
of ways to do this is: 



(N — X m )\ ' 



6 



(ii) Now we take the Y m Y-mers in the (m — l)-th layer 
that are extending down to the lower layer. Number of 
ways of extending these to the (N—X m ) unoccupied sites 
below is clearly 

(N - X m )\ 



(N - X n 



Y m )\ 



(iii) Connect the remaining N — X m Y-bonds between 
layers m — 1 and m to sites in layer m not yet connected 
by X bonds. The number of ways to do this is 

(N — X m )\. 

(iv) Repeat the last procedure with the N — Y m re- 
maining Y-bonds between the layers m and m + 1. The 
number of ways to do this is 

(N - Y m )\. 

(v) Finally, the (N — X m — Y m ) sites in layer m, that 
are unoccupied so far, are divided into three groups: x m 
topmost sites of new X-mers, y m topmost sites of new 
Y-mers, and the unoccupied sites. Clearly the number of 
ways to do this is: 

(N-X m -Y m )\ 



x m } -y m \(N - X m - Y„, 



VrnV- 



The product of these factors gives the total number of 
ways of adding the m-th layer as: 



N\(N-X m )\(N-Y m )\ 



x m \y m \(N - X n 



y m V- 



(14) 



Finally, multiplying these factors for different m, we get 
the configurations with specified {x m ,y m } have total 
number of configurations C({x m ,y m }) given by 



C({x m ,y m }) = Y[ 



N[(N-X m )l(N-Y m )\ 



, %rn !y?n ! (-A^ X m Yjn Vrnj- 
rn—1 a v ' 

Putting in the corresponding activity factors, the grand 
partition function for the whole lattice with M layers is 



(N\)2M-2 
V '' {lm,h} 



C({x m , y m }) z i 



where the sum is over all possible values of {x m ,y m }), 
and we have divided the multiplicity factor by the num- 
ber of configurations of the random lattice (7V!) 2M ~ 2 , to 
get the average partition function. Note that Z av = 1 for 
Z\ = Z2 = 0, as expected. 

The summation over {x m ,y m } yields at most a factor 
of order N 2M . Since the summand is of order exp(iVM), 
for large TV, we can ignore the summation over {x m , y m }, 
and replace the summation with the largest term, with 
negligible error. For the summand to be maximum with 
respect to Xj, we set 



y m })zi 

C({x 

in : Dm}) 



1. 



(17) 



This gives 
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s=0 



(N - X j+S - x j+s - Y J+S - y j+s ) _ xj + 1 



(N - X j+S ) 



and similarly 
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(N - X j+S - x j+s - Y j+a - y j+s ) _ (yj + 1) 



(N-Y j+S ) 
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Writing the maximising values as f* — Xj/N and g* = 
yj/N, the equations satisfied by f* and g* in the limit 
N tends to infinity are 
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•' J 



s=0 
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p X (j + s) + fj+s l-p*U)+fj 
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z TT 1 ~ Pxjj + S) - py(j + s) = 

Z2 Lo 1 - Pv ( j + s ) + 9*i+s 1-Pyti)+ 9* 



(21) 



where p x (rn) and p v (m) are the fractions of sites in layer 
m covered by X-mers and Y-mers respectively. Clearly, 



fc-i 



fc-i 



Px (J) = Y1 fj-*> Pv (■?) = X] ■' 



(22) 



s=0 



s=0 



These equations connect /* and g* to their value of 
fj+ s and g* +s , with s = 1 to A; — 1. These may be consid- 
ered as recursion equations for /*,<?*• These recursions 
work in the direction of decreasing j. 

These equations have a simple interpretation. In the 
equilibrium state, Z\ is the ratio of the probability that 
a randomly chosen site will be the head of a X-mer to 
the probability that an X-mer can be placed with this 
site as the head. The probability that the chosen site is 
empty is [1 — p x (j) — Py{j)\- Given that a given site is 
empty in layer (j' — 1), the conditional probability that 
the site connected to it in the layer j' by an X-bond is also 
empty is ^-p^j^/'* ' ■ Multiplying these probabilities 

for (k — 1) consecutive layers, we get the probability that 
a given site in layer j can be the head of an X-mer to be 



ntoti 



Px{j + s) - Py{j + s)\ 



ntl^ -px(j + s) +/?+.)] 

The probability that the chosen site in the jth layer is 
the head of an X-mer is /*. The ratio of these is 1/zi, 
which gives Eq. ([T^]). 

The simplest solution of this is a fixed point solution 
with f* — /*, g* — g* , independent of j (away from the 
boundaries). Then f*,g* satisfy the equations 



zi(l 
zi{\ 



kf* 
kf 



kg 
kg 



k 



f [1 



(k 



1)/' 
1)5 



ife-i 



*ifc-i 



(23a) 
(23b) 
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These equations are the same as Eq. (fTU|) , and we have 
recovered the Bethe approximation. However, for the 
RLTL lattice, the limit of fully packed lattice is well- 
defined, and causes no difficulties. 

From Eq. (fl"5]) . the entropy per site (divided by ks) is 
easily seen to be 



, k- 1 \, / k- 1 

S{Px, Py) = [ 1 — Px I In I 1 — p x 



fe- 1 \, / fc- 1 

1 - — — Py In 1 - — — p y 



-(l-p)ln(l-p)-^ln^--^ln^(24) 

where p = p x + p y is the total density. The same ex- 
pression for entropy was obtained by DiMarzio [2l|, |4(| 
who used an approximate counting technique for count- 
ing configurations on cubic lattices in any dimension. 
Also, Eq. (j2~4")l coincides with the expression for entropy 
that one obtains by using Gujrati's prescription for cal- 
culating free energies on the Bethe lattice [37| . 

It is easy to see that this expression for the entropy 
is not everywhere convex. When the value of s(p x ,p y ), 
calculated as above, turns out to be in the non-convex 
region, it is easily seen that a much larger contribution 
to the partition function comes from {x m ,y m } that are 
not nearly uniform, and in a canonical ensemble at fixed 
p x and p y , the lattice will show phase separation, with 
one region having higher density than the other. The net 
effect of this to replace non-convex parts of the entropy 
function by a convex envelope construction. Thus, we 
write the true entropy s(p x ,p y ) as 

s(Px,P y ) = CS [s(p x ,p y )\ , (25) 
where C£ denotes convex envelope. 



V. ISOTROPIC-NEMATIC TRANSITION 

In this section, we analyse the isotropic-nematic tran- 
sition when the coordination number is four. From 
Eq. (|24p , the expression for entropy, it is straight forward 
to determine the ordering that has maximum entropy. 

Consider the system at a fixed density p = p x + p y . 
We define the order parameter ip by 



Px - Py 
P 



(26) 



Then, it is easy to study the variation of s(p x ,p y ) as a 
function of ip, for fixed p. We find that for small p, the 
entropy has a single maximum at ip = 0, but for large 
enough p > p c , it develops two symmetrically placed 
maxima (see Fig. [5]). For small ip>, we can expand the 
entropy in a power series around ip = as in the stan- 
dard Landau treatment 




0.36 1- P < Pc 

0.34 
0.32 
0.3 
0.28 
0.26 
0.24 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
V 

FIG. 5. Entropy a function of the order parameter ip 

for different densities. The data are for q — 4, k = 4 when 
the transition is continuous. Entropy has one peak for small 
densities and two symmetric peaks for large densities. The 
dotted lines denotes the convex envelope. 



The coefficient of the quadratic term B(p) changes sign 
at p = pc, and is negative for larger p. Thus, for p > p c , 
we have a nematic phase, with nonzero value of ip. The 
critical exponent for (3 for the order parameter takes the 
classical Landau theory value 1/2. 

The value of p c can be determined easily from Eq. (|23[) . 
For z\ = Z2 — z, we note that both /* and g* are solu- 
tions to an equation of the form 



k-1 , 
1 — x — constant. 



(28) 



The left hand side of the equation is a function of x 
that starts at for x = 0, increases to a maximum value, 
and then decreases monotonically, and reaches a positive 
finite value (l/fc) fc_1 , for x = 1. When the right hand 
constant is small enough, there is only one real valued so- 
lution of this equation, and /* = g* . For a range of values 
of the constant, there are exactly two distinct solutions. 
At the critical point, the two solutions are degenerate. 
This occurs where the function is maximum, i.e., at 



1 



1 



(29) 



Then at this point /* = g* = x*/k, and p x — p y — 
l/(k — 1). Correspondingly, we have p c = 2/(k — 1). The 
corresponding value of critical activity, from Eq. (|23[) . is 



(k-1) 



2fe-2 



[k(k - 3)] fc 



(30) 



s(p x ,p v ) = A(p)-i; 2 B(p)+i; i C(p) + 



(27) 



We note that the value of z c is finite only for k > 4. In 
Fig. |6l we show the variation of the modulus of the order 
parameter with density p for different rod lengths 
k. | ip\ is non-zero for densities larger than the critical 
density. 



8 




0.1816 



0.1812 - 



0.1808 




FIG. 6. The modulus of nematic order parameter ip as a 
function of the density p for different rod lengths k. The data 
are for q = 4. 



VI. LATTICES WITH COORDINATION 
NUMBER q > 6 

The analysis of q > 6 is very similar. It is easy to 
check that for general even q, the fixed point solution 
which is independent of the layer index j, satisfies the 
self-consistent equations 



kz(l - pf = Pi 



1 - 



k - 1 



k 



-Pi 



fc-i 



i=l,. 



,f, (3D 



where i labels the q/2 different bond types, and p = pi 
is the total density of sites that are part of fc-mers. The 
entropy per site generalises to 



q/2 



_(l-p)hi(l-p)-X;?ln?. 



(32) 



The low density phase is isotropic, with same for 
different i. We define the order parameter to be ip = 
(pi — P"i)l P% with p 2 = ... = Pq/2- But now the entropy 
function has no symmetry under ip — t —tj). Then, the 
expansion of s(pi) in powers of ip contains cubic terms. 
Fig. [7] shows the behaviour of entropy as a function of ip 
for different densities p when q = 6. For small p, there 
is a single maximum at ip equal to zero. For larger p, a 
second local maximum at a non-zero ip appears, and at 
some value of p this becomes of equal height. Then the 
order parameter jumps discontinuously, as the density is 
increased. 

In Fig. [51 we show the variation of the order parameter 
ip with density p for q = 6 and different values of k, and 
for k = 5 and different values of q. The first order tran- 
sition is clearly seen, with the critical activity increasing 
with q and decreasing with k. 



FIG. 7. Entropy a function of the order parameter ip 

for different densities. The data are for q = 6, k = 8, when 
the isotropic-nematic transition is first order. The different 
curves correspond to (a) p = 0.3890, (b) p = 0.3910 w p c , (c) 
p — 0.3930. The dotted line shows the convex envelope. 
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0-6 w1 ;q 1 

^_ k-5;q-6 
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0.2 0.4 0.6 0.8 1 
P 



q=6; k=4 
q=6; k=5 
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P 

FIG. 8. The order parameter as a function of density p for 
q = 6 and varying k. There is a first order transition at a 
critical value p c which decreases with k. Inset: The order 
parameter ip as a function of p for k — 5 and varying q. p c 
increases with q. 



For coordination number greater than four, it is not 
possible to determine the critical density exactly. The 
critical densities obtained by numerically comparing the 
entropy of the isotropic and nematic phases are sum- 
marised in Fig.[S] For a fixed value of q, p c ~ hi(q)/k for 
large k. 

As expected in discontinuous transitions, in a range 
of values for the densities around the critical density p c , 
entropy has a local maximum for both the isotropic (ip = 
0) and nematic phases. In Fig.[TUl the values of densities 
at which these local maxima appear and disappear, along 
with the critical density are shown for k — 5 and k = 8 for 
different values of q. Only for q = 4, where the transitions 
are continuous, all the densities coincide. 



9 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 





q=4 


+ 


G V\\ 


q=6 


X 




q=8 






k q= 10 


□ 




.Xv q=12 


■ 




xx\x q= 14 


o 




q=16 


• 




\\"$%s^ q=18 


A 









10 



FIG. 9. The critical density p c of the isotropic-nematic phase 
transition as a function of k for different values of the coordi- 
nation number q. For 5 = 4 the transition is continuous. 




FIG. 10. Coexistence and limits of stability values of density 
of the isotropic and nematic phases as functions of q for two 
different rod lengths. For given values of q and k, the lowest 
point is the is the density at which s(tp) develops a second 
local maximum, the intermediate point is the critical density 
p c , and the highest point is the density at which s(ip) is no 
longer a local maximum at tp — 0. The lines are guides to the 
eye. 



It is possible to determine the spinodal density p s , the 
density at which the entropy s(i/j) no longer has a maxi- 
mum at if) = 0. As in the analysis of q = 4, the spinodal 
is still given by the condition x* = l/(k — 1). But now 
that there are q/2 distinct directions, the density at the 
spinodal point is q/[2(k — 1)]. If q < 2(/c — 1), the spinodal 
density is less than one, and entropy will have a maxi- 
mum at the nematic fixed point. Thus, there is always a 
phase transition for q < 2(k — 1). 

We now address the question of whether for coordina- 
tion number q, there is a minimum rod length k m i n below 
which the isotropic-nematic transition is absent. From 



TABLE I. The minimum rod length k m i n required for an 
isotropic-nematic transition as a function of coordination 
number q. 
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the analysis of the spinodal, 



k n 



< 



(33) 



The precise value of k m i n is computed numerically and 
the results are shown in Table |U We find that k m i n ~ 
In(g). 



VII. DISCUSSION 

We are not able to address the nature of the second 
transition in this study, as on the RLTL lattice, there is 
none. However, on this lattice, the limit of fully packed 
lattice is quite interesting. We find that in this limit, the 
system has long-range nematic order, but the ordering is 
not complete, and there are small islands of "wrongly" 
oriented fc-mers in a sea of aligned fc-mers. The small 
concentration of these wrongly oriented rods is entropi- 
cally stabilised. 

The order parameter tp is easy to determine for k = 4, 
when 



^ w 0.96225, 
3(v f~ 4) « 0.93031, 
±f « 0.86667, 



k=4, q=4, 
k=4, q=6, 
k=4, q=8. 



(34) 



The limiting value of ip grows monotonically with k, be- 
ing equal to 1 in the limit k — > oo. It is easy to see 
that the nematic order parameter has the asymptotic be- 
haviour 



1 -ip 



<l 



2k 



k-l 



k ->• 



(35) 



We can also look for a periodic solution of period k, 



where x m+k = x m ,y m+k 



In this case, the 2k 



independent parameters fg,gti with s 
the equations 



to k— 1 satisfy 



where 



- pi - 

- 92 

Pi - 



fe-i 

s=0 



zi{l 

22(1 



P2 



- Pi 

- Pi 

k-l 
s=0 



P2TIA, 

P2) k /B, 



(36) 



(37) 
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and 



fc-i 



k-l 



A = I] I 1 - Pi + f* ] ; B = II t 1 - + tfl ■ (38) 



In Eq. (|37|) . the left hand side is a function of the form 
(x + cir 2 ), with c positive, and the right hand side is in- 
dependent of j. Hence, the only positive real solution of 
this equation is of the form f* independent of j, and we 
do not get a non-trivial periodic solution. Note that a pe- 
riodic solution would correspond to smectic-like layered 
ordering, and our solution rules it out. 

We can easily extend our treatment to semi-flexible 
rods, where all the rods are aligned in the direction of 
increasing layer number, but a fc-mer lying on an X-bond 
between layers j and (j = 1) can bend, and lie on a Y- 
bond between layers + and (j+2), with some energy 
cost. Solutions for flexible and semi-flcxiblc polymers on 



Bethe and Husimi lattices can be found in Refs. [41|, 42 1. 

There are other models like the ANNNI model 43J, 
where exact solution in dimensions greater than one is 
not possible, and the equilibrium state shows spatial 
structure. These are usually discussed in the spatially 
varying mean-field approximation. We feel that study- 
ing such model on RLTL lattice can take into account 
the short-range correlations in these systems better. 
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